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^ ' A new rigorous light-cone path integral approach to the Landau-Pomeranchuk-Migdal 

•/^ I effect in QED and QCD is discussed. The rate of photon (gluon) radiation by an electron 
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(quark) in a medium is expressed through the Green's function of a two-dimensional 
Schrodinger equation with an imaginary potential. In QED this potential is proportional 
to the dipole cross section for scattering of e^e" pair off an atom, in QCD it is proportional 



Oh! 

^ • to the cross section of interaction of the color singlet quark-antiquark-gluon system with 



a medium constituent. In QED our predictions agree well with the photon spectrum 
measured recently at SLAG for 25 GeV electrons. In QGD for a sufficiently energetic 
quark produced inside a medium we predict the radiative energy loss AEg oc L^, where L 
is the distance passed by the quark in the medium. It has a weak dependence on the initial 
quark energy Eg. The dependence transforms into as the quark energy decreases. 
We also give new formulas for nuclear shadowing in hard reactions. 



E-mail: bgz@landau.ac.ru 



1 Introduction 



In 1953 Landau and Pomeranchuk [jT| predicted within classical electrodynamics that 
multiple scattering can considerably suppress bremsstrahlung of high energy charged par- 
ticles in a medium. In the high energy limit they obtained the photon radiation rate 
oc l/^/tJ {uj is the photon frequency), which differs drastically from the spectrum for an 
isolated atom cc l/u. Later, this result was confirmed by Migdal 0, who developed a 
quantum-mechanical theory of this phenomenon. The physical mechanism behind the 
suppression of the radiation rate in a medium is the loss of coherence for photon emission 
from different parts of the charged particle trajectory at the scale of the photon formation 
length. 

Since the studies by Landau and Pomeranchuk and Migdal 0, the suppression of 
radiation processes in medium, called in the current literature the Landau-Pomeranchuk- 
Migdal (LPM) effect, has been studied in many theoretical papers [3-21]. However, only 
recently the first quantitative measurement of the LPM effect for high energy electrons was 
performed at SLAG [^. Altogether, this experiment corroborated LPM suppression of 
bremsstrahlung. Unfortunately, contamination of the SLAG data by multiphoton emission 
makes it difficult to perform an accurate comparison with the theoretical predictions in 
the entire range of photon energies. Nonetheless, the experimental spectrum for a thin 
gold target at photon energies from 5 MeV up to 500 MeV for 25 GeV electron beam, 
where the multiphoton emission gives a small contribution, agrees with predictions of Ref. 



Taj. The calculations of Ref. |15] were carried out within a new light-cone path integral 



approach to the LPM effect which we developed in Ref. [14|. In Ref. [1^ the radiation 



rate was calculated including inelastic processes and treating rigorously the Goulomb 
effects. In all previous analyses the inelastic processes were neglected, and the Goulomb 
effects were treated in the leading-log approximation. This approximation works well in 
the limit of strong LPM suppression in an infinite medium. However, it is not good for 
real situations because of the uncertainty in the value of the Goulomb logarithm. 

The approach of Ref. is also applicable in QGD. Analysis of the LPM effect 



in QCD is of great importance for understanding the longitudinal energy flow in soft 
and hard hadron-nucleus collisions and the energy loss of a fast quark produced in deep 
inelastic scattering on nuclear target. It becomes especially interesting in connection with 
the forthcoming experiments on high energy AA-collisions at RHIC and LHC, where the 
formation of quark-gluon plasma (QGP) is expected. The energy loss of high-p^ jets 
produced at the initial stage of AA-coUision may be an important potential probe for 
formation of QGP. The first attempt to estimate the radiative quark energy loss AEq in 
QGP was made by Gyulassy and Wang |T^. They modelled QGP by a system of static 
scattering centers described by the Debye screened Coulomb potential oc exp{—rfiD)/r, 
where /i_D is the color screening mass. The authors studied emission of soft gluons in 
the region of small transverse momenta fc^ << fiD, which, however, gives a negligible 
contribution to AEq. Analysis of soft gluon radiation without a restriction on the gluon 
transverse momentum in the limit of strong LPM suppression within the Gyulassy- Wang 



(GW) model was performed in Refs. |TT|, A similar analysis for cold nuclear matter 



was given in Ref. ^7j. However, the authors of Refs. |TT|, ^ used some unjustified 
approximations. For instance, the quark-gluon system emerging after gluon emission was 
treated as a pointlike color triplet object. Q A rigorous quantum treatment of the induced 
gluon radiation was given for the first time in Ref. [0 (see also 

In the present paper we discuss the LPM effect in QED and QCD within the approach 
of Ref. |T^. We give special attention to technical details omitted in our previous short 
publications. We also consider from the viewpoint of LPM suppression nuclear shadow- 
ing in hard reactions. New formulas for shadowing that take into account the parton 
transverse motion are derived. 

The presentation is organized as follows. Section 2 is devoted to the LPM effect in 
QED. Using the unitarity of scattering matrix we express the cross section of photon emis- 



^ Note that, after submission of the present paper in November 1997 to Phys. At. NucL, R.Baier, 
Yu.L.Dokshitzer, A.H.Mueller and D.Schiff (iiep-pii/9804212) had reanalyzed the induced gluon radiation 
without using the approximation of the pointlike qg system. 
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sion through radiative correction to the transverse electron propagator. It is calculated 
within time-ordered perturbation theory (PT) in coordinate representation. The radiation 
rate is expressed through the Green's function of a two-dimensional Schrodinger equation 
with an imaginary potential which is proportional to the dipole cross section for scattering 
of e~'"e~7 state off an atom. We demonstrate that the cross section of photon emission 
can be written in the form analogous to the Glauber amplitude for elastic hadron-nucleus 
scattering. This representation allows one to view LPM suppression as an absorption 
effect for e"'"e~7 system. We compare the theoretical predictions with the data of the 
SLAG experiment In section 3 we discuss the LPM effect in QGD. As in QED, the 
radiation rate is expressed through the Green's function of a two-dimensional Schrodinger 
equation. The corresponding imaginary potential is proportional to the total cross section 
for a three-body quark-antiquark-gluon state. We evaluate the quark energy loss for cold 
nuclear matter and QGP using the oscillator parametrization for the imaginary potential. 
For a high energy quark incident on a nucleus we predict AE'^ ~ 0.1-E'q(L/10 fm). For a 
fast quark produced inside a medium we obtain /S.Eq oc (L is the quark path length 
in the medium) while at sufficiently small energy AE'g oc L. In section 4 we discuss the 
LPM effect for hard reactions on nuclear targets. We demonstrate that our approach to 
the LPM effect can be used for an accurate evaluation of nuclear shadowing. In section 5 
we summarize our results. 

2 The LPM effect in QED 

2.1 General expression for the radiation rate 

We begin with the LPM effect for bremsstrahlung of a fast electron. We consider an 
electron incident on an amorphous target of a finite thickness. Multiphoton emission will 
be neglected. The probability of photon emission, P^, is connected with the probability 
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Pe to detect in the final state one electron by the unitarity relation: 

Pe + P^ = l. (1) 

In the absence of interaction of the electron with the quantum photon field, we have 
Pe = I. Consequently, Eq. (|lD can be rewritten as 

P, = - {6Pe - 5Pr) , (2) 

where 5Pe is the radiative correction of order a [a = 1/137) to Pe- On the right-hand side 
of Eq. d^) we subtracted the vacuum term, which takes into account the renormalization 
of the electron wave function for initial and final electron states. 

We will evaluate 5Pe in time-ordered PT P3| , The corresponding matrix element 
is generated by transitions e — e'7 e. In a medium the electron transverse momentum 
is not conserved. For this reason it is convenient to use the coordinate representation of 
time-ordered PT in which the transitions e — > e'7 e will reveal themselves through the 
radiative correction to the electron wave function. Let us first consider the wave function 
of a fast electron neglecting interaction with the quantum photon field. In the vacuum 
the radiative-correction-free wave function of a relativistic electron with longitudinal mo- 
mentum Pz » me {me is the electron mass) can be written as 

ipit, r) = exp[-ipz{t - z)](f){t, p) , (3) 

where r = {z,p), and the time-dependence of the transverse wave function 4>(t,p) is 
governed by the two-dimensional Schrodinger equation 

with the Hamiltonian 



i70(t,p), (4) 



H=^^±^. (5) 

Here q is the operator of transverse momentum, and the Schrodinger mass is /ie = Pz- 
Eqs. P) hold for each helicity state. At high energy the electron propagates nearly 



along the light-cone t — z = const, and in Eq. (^) the variable t can be viewed as the 
longitudinal coordinate z. For this reason, we will henceforth regard the transverse wave 
function <^ as a function of z and p. Eq. (^) allows one to write the following relation, 
connecting (j){z, p) at planes z = Zi and z = Z2, 

(P{Z2, P2) = J dplKe{p2, Z2\pl, Zi)(f){zi, p^) , (6) 

where 

Hiie{p2-Pif im^Az] 



ire(P2,-2|Pi,.i) = ^exp 



2Az 2/ie 

is the Green's function of the two-dimensional Hamiltonian (j^), with Az = Z2 — Z\. At 
high energies spin effects in interaction of an electron with an atom vanish, and equations 
analogous to (|^) and (|^) hold for propagation of an electron through the medium as well. 
The corresponding propagator can be written in the Feynman path integral form 

where p = dp/dz, and U{p, z) is the potential of the medium. 

The photon wave function can also be written in the form (^. Using the representation 
(0) for the electron and photon wave functions, we can obtain for the radiative correction 
to the transverse electron propagator associated with e'7 intermediate state 



Ke{p2, Z2\Pi, zi) = J Vpexpli J dz 



Z2 



5Ke{p2,Z2\pi,Zi) = - J dx J d^i J d^2 J dTidT2gi^u^2,x) 

zi 5i 

'XKe{p2, Z2\t2, 6)^e'('r2, 6 1 ^1 , 6)^7 (^^2, 6 1 ^1 , 6)^e (Ti , ^1 1 Pi , Zi) . (9) 

Here the indices e' and 7 label the electron and photon propagators for the intermediate 
e'7 state. The Schrodinger masses that appear in the Green's functions Ke' and 
are /Xg' = (1 — x)fie and /i^ = xfie, where x is the light-cone fractional momentum of 
the photon. The vertex operator g{C,i,^2,x), including all spin effects associated with 
transitions e —>■ e'7 —>■ e, is given by 



a [4 — 4x -|- 2x^1 ,^ , ,^ , amlx 



^(6, 6, x) = ^ ^ v(6) • v(ei) + , (10) 
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where 

and Ve' are the transverse velocity operators, which act on the corresponding propaga- 
tors in Eq. (P). Two terms on the right-hand side of Eq. ( p!OD correspond to the e — > e'7 
transitions conserving and changing the electron helicity. 

Now we have all ingredients necessary for calculating the radiation rate. We consider 
a target with a density n{z) independent of the impact parameter, and assume that n{z) 
vanishes as 1^1 — > 00. In terms of the transverse Green's function Kg and 6 Kg the radiative 
correction 6Pe can be written as 

5Pe = 2Re J dp^dp\dp2(l){zi, p^)(p*{zi, p[){5Ke{p2, Z2\pi, zi)Kl{p2, Z2\p'i, zi)) , (11) 

where (...) means averaging over the states of the target, and the points Zi and Z2 are 
assumed to be at large distances before and after the target, respectively. The initial 
electron wave function is normalized by the condition 

'dp\<p{zup)\' = l . (12) 

For a high energy electron one can neglect the longitudinal momentum transfer associ- 
ated with the interaction with a medium potential. For this reason the unitarity relation 
(0) is also valid in the differential form in the light-cone variable x. Then, using Eqs. (^, 
(y) and (p!TD, we can obtain for the radiation rate [we suppress the vacuum term, which 
will be recovered in the final formula (|33D ] 
dP 

= 2Re J d^i J d^2 J dpidp[dTidT[dT2dT2dp2(f){zi, pi)(f)*{zi, p[) 

X5(6, 6, X)S{P2, p2, Z2\t2, T2, 6)^(t2, 61^1, t[,^i)S{Ti, t[, ^i\p^, p[, Zi) , (13) 

where 

5(P2,P2,6|P1,P;,6) = (ire(P2,6|Pl,6)^:(p'2,6|p'l,6)) (14) 

is the evolution operator for the electron density matrix in the absence of interaction with 
the photon field, and 

M(p2,P2>6|Pl,P'l,6) = (ire'(P2, 6)^7(^2, 6|Pl,ei)^:(P2,6|p'l, 6))- (15) 



In deriving Eq. (|T3|) we used the convolution relation 

Ke{p2,Z2\pi,Zi) = J dTidT2Ke{p2, Z2\T2,^2)Ke{T2,^2\ri,^l)Ke{ri,^l\pl, Zi) . (16) 

Using the path integral representation for the transverse Green's functions, we can 
rewrite Eqs. ([T^), ([T5|) in the form 



H{Pe}.{p'e}) 



(17) 



M{p2,p2,^2\pi,p'i,^i) = J Vp^,Vp^Vp^exp 1^ J d^ifie'pl + ^i^pt^ - Hepl) 

^(6-6; 



Here we introduced the photon formation length 



Lf = 2/2^ 



1 — X 



2^e(l - x) 



mix 



<f({Pe'},{Pe}) ■ (18) 



(19) 



which appears in Eq. (|TBp owing to the difference between the phases of the wave functions 
of e and e' states. Notice that the value of L/, emerging in deriving Eq. (|15), agrees 
with the estimate based on the uncertainty relation AEAt ~ 1. In the following we 
will assume that Lf is much larger than the atomic size. The boundary conditions for 
trajectories in Eq. are pMi,2) = Pi,2, Pe(^i,2) = p'i,2) and in Eq. (|18D Pe/,^(6,2) = Pi,2, 
Pe(^i,2) = p'l 2- The phase factor <l> in Eqs. ( p^Tf ) and ([T8|) that takes into account 
interaction with a medium potential is given by 



^({PJ,{P,}) = (expj^e J daU{p,{0,0-Uip,iO,0] 



(20) 



Note that this phase factor can be viewed as the one for propagation through the medium 
of e+e" system. However, it should be borne in mind that the "positron" kinetic energy 
term in Eqs. ([T7|), ([T^) is negative. 

We will neglect the correlations in the positions of medium atoms. In this case 
$({Pj}, {Pj}) can be written as 

1 



<f({p.},{p,}) = ^-^ /«0 



X 



I rfb^l-expjze I deW.(O-b,e'-0-^(P,(O-b,e'-0] 



N 



(21) 



where (p{r) is the atomic potential, is the number of the atoms in the targets, ( ... )a 
denotes averaging over the states of the atom. After exponentiating Eq. (|21|) can be 
written in the form 



H{P.}, {P,}) = exp [-^ / d^n{Oa{\p,{0 - P,(OI) 



(22) 



where 



a(\p\) =2jdh(l- exp jze J d^[^{p - b, - V^(b, 0]})^ (23) 

is the dipole cross section for scattering of e^e~ pair of the transverse size p on the atom. 
In arriving at (|22| ) we neglected the variation of |Pj(0 ~ PjiOl ^^e longitudinal scale 
on the order of the atomic radius a. This is a good approximation for Lf » a. 

For an atomic potential (f{r) = (Ze/47rr) exp(— r/a) (a ~ tqZ^^^^, is the Bohr 
radius) a{p) in the Born approximation is given by 



a \a 



(24) 



where Ki is the Bessel function. For p <^ a, which will be important in the problem under 
consideration, we have a{p) ~ C{p)p'^ , where 



C{p) = A7r{Zaf 



\ (2a\ (1-27) 



7 = 0.577. 



(25) 



For nuclei of finite radius Eq- (0) holds for p > Ra, and C{p < Ra) = C{Ra)- In 
section 2.4 we will give a more accurate formula for C{p), which will be used in numerical 
calculations. 

The phase factor (p^) is independent of (pj + Pj)/2. This allows one to calculate the 
path integral (0) analytically. This gives the following result p5| 

J ipe 



'5'(p2,p2,6|Pl,pi,^l 



Pe 



exp 



I \2l 



[(Pl - P2)' - (P'l - P2) 



(26) 



TsiO = (Pi - P'i)%# + (P2 - P2)^^TT^ , = 6 - 6 • 



Ae 
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In the case of the integral (|I8|) we introduce the Jacobi variables a = (/ie'Pe/+/i-yP^)/(/ie' + 
/i^) and p = Pg, — . Then after analytical path integration over cx{^) and Pe(0 Eq. 
fllSl) we arrive at 



M{p2,p2,^2\Pl,Pl,^l 



X j Vpexp Y J 



(Pi - P2) 



IPl - P2) 



Traii) = (Pl - Pi) + (p2 " P2) + 



2 



(27) 



Ae ' ' (/ie'+Zi^)' 

where fie'-y = /^e'/^7/(/^e' + ^^7) = Eex{l — x) is the reduced Schrodinger mass of the e'7 
system. It follows from Eqs. (0), (0) that for the factors S and M there hold relations 



/ 0?p2'5(p2,p2,6|pi,Pi,6) ='5(Pi-pi), 
I (ip2M(p2,P2,6|pi,p'i,ei) =5(Pi-p'i)/C(0,6|0,e 



where 



^(p2.6|Pi,6) = J Vpexpli J 



/ie'^P^ , ,^(0^(|P|^) 



+ z- 



(28) 
(29) 

(30) 



2 2 

is the Green's function of a two-dimensional Schrodinger equation with the Hamiltonian 

.2 



n 



q 



2/i, 



e'7 



^^(P,0 



+ ^(p,0, 
.^(O^(lpk) 



(31) 



(32) 



The Hamiltonian (0) describes the electron-photon relative transverse motion in the e'7 
system. The form of the imaginary potential ( P^ reflects the fact that after integrating 
over p2 in Eq. ( pOj) the "positron" trajectory Pe(0 coincides with the trajectory of the 
center-of-mass of the e'7 system. 

Substituting (p6|), (|2^) into (|T3|) and integrating in (|1^) over the transverse variables 
with the help of (^), ([29|), along with the normalization condition (|12|), we finally obtain 
(we set —Zi = Z2 = 00, and recover the vacuum term) 

^ = 2Re 1 ^6 / f/6exp(-^j(7(ei,e2,a;)[/C(0,6|0,ei)-/C,(0,6|0,ei)] • (33) 
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Here JC^ is the Green's function for the Hamiltonian (pl|) with v{p,^) = 0. 

Thus, we expressed the intensity of the photon radiation through the Green's function 
of the Schrodinger equation with the imaginary potential (|3^) . Notice that the dipole 
cross section in Eq. (^) can be viewed as an imaginary part of the forward scattering 
amphtude for the three-body e"^e~7 system. In our analysis we neglected interaction 
of the photon with atomic electrons, which becomes important only for extremely soft 
photons [|], 0. The inclusion of this interaction leads to appearance of a real part in the 
potential 

It is worth noting that an equation analogous to ( p3D holds also in more general case 
of photon emission on a random external potential if the photon formation length is much 
larger than the potential correlation radius. In this case the integral over ^ in Eq. (pO|) 
defines a Gaussian random quantity. Using the formula {exp{iA)) = exp(— (A^)/2), which 
is valid for a Gaussian random quantity, we find that the phase factor (pOD takes the form 



(34) 



V{p, = -2ma J d^'{AU{p, OU{p, O) , (35) 

— oo 

where AU (p, — U{p,^) — U (0, The function (^) will play the role of the imaginary 
potential in the Hamiltonian (|3lD. 



2.2 Bremsstrahlung in an infinite medium in the oscillator ap- 
proximation 

To proceed with analytical evaluation of the radiation rate we take advantage of the 
slow p-dependence of C{px) at px < l/mg, which, as will be seen below, are important 
in Eq. (|33D. Evidently, to a logarithmic accuracy we can replace (|32D by the harmonic 
oscillator potential with the frequency 
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Here pe// is the typical value of p for trajectories dominating the radiation rate. Making 
use of the oscillator Green's function 



Kosc{P2, Z2\Pl, Zi) 



■ exp 



if^iipl + pI) cos {n^z) - 2P1P2] 



(37) 



2m sin {Vt^z) \ 2 sin {Vt/^z) J 

after some algebra one can obtain from Eq. (|33|) the intensity of bremsstrahlung per unit 
length in an infinite medium 



dxdL 



n 



X] 



\C{1/ 



rrir-. 



'da 
dx 



BH 



'da 
dx 



BH 



(3^ 



nf \ / sf 

where 1] = Lf\Q\. In Eq. ( ^Sf ) we factored out the Bethe-Heitler cross sections conserving 
(nf) and changing (sf ) the electron helicity, which to a logarithmic accuracy can be written 
as 



'da 
dx 



BH 



aC{l/me) [4-4x + 2x2 



nf 



'da' 
dx 



BH 



X 



aC{l/me 



sf 



■X . 



(39) 
(40) 



The factors Snf, Ssf in Eq. (|38| ) are given by 



Snfiv) 



r]V2 



y 



Ssfiv) 



6 



yj \ r]V2, 
1 1 \ 



cos 



+ sin 



dy \ \ — I exp , _ 

r ^ \y shy J \ r]V2 



y \ . y 

sm 



(41) 



(42) 



At small rj Snfirj) ~ 1 — 16ri'^/21, Ssf{i]) ~ 1 — 317^^/21, and the Bethe-Heitler regime 
obtains. Up to the factor C{peffx)/C{l/me), which is slowly dependent on t], the sup- 
pression of bremsstrahlung at ^ 1 is controlled by the asymptotic behavior of the 



suppression factors (^TD 



(43) 



The value of pe// can be estimated with the help of Eqs. (^l)) , (^21) . The variable of inte 



gration in (^T]), (42) in terms of in Eq. ( p3l) equals |A^fi|. Therefore, for typical value 



of contributing to the integral (0), ^C,eff, we have A,^, 



eff 



min(L/,l/|fi|). 
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Note that L'j plays the role of the effective medium-modified photon formation length. 
Having ^ieff we can estimate Pe// from the obvious Schrodinger diffusion relation: 

Pe//~(2Aee////ie',)'/'. (44) 

In the low-density limit, when ^ 0, this relation yields pe// ~ l/iTieX, and the right- 
hand side of Eq. (|38D goes over into the Bethe-Heitler cross section times the target 
density. In the soft photon limit (x —>■ 0) at fixed n, rj becomes much greater than unity. 
In this regime of strong LPM suppression, using the asymptotic formula for Snj, one can 
obtain from Eqs. (^8|), ( ^9l) 



dxdL V nEeX 

r 1^1/4 
with peff ~ n{ZaynEeX^ log(2/aZi/3) . This resuh agrees with Migdal's prediction 

obtained within the Fokker-Planck approximation in the momentum representation. 

Our suppression factors (^), ( ^21) also agree with those obtained in Ref. 0. Equivalence 

of the the oscillator approximation in coordinate representation to the Fokker-Planck one 

in momentum representation is not surprising. Making use of Eq. (|26|) one can easily 

show that cr(p) oc p^ leads to the Gaussian diffusion in the momentum space. That is, we 

have a diffusion described by the Fokker-Planck equation. 

The oscillator approximation simplifies greatly evaluation of the radiation rate. It 

allows one to obtain simple formulas for suppression factors for a finite-size target as 

well. The corresponding analysis within Migdal's approach was performed in Ref. 

Unfortunately, the oscillator approximation is accurate only for strong LPM suppression, 

when one can neglect the variation of the factor C(p). This effect must be taken into 

account to evaluate accurately the radiation rate in an infinite medium in the regime 

of small LPM suppression and for finite-size targets. In the next section we represent 

Eq. ( p3D in a different form which is more convenient for numerical calculations with a 

rigorous treatment of the Coulomb effects. 
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2.3 Glauber form of the radiation rate 



In this section we demonstrate tliat Eq. ( |33D can be rewritten in a form analogous to 
the Glauber amplitude for elastic hadron-nucleus scattering. Let us expand the Green's 
function /C in Eq. (RSl) in a series in the potential v 



^(P2) Z2\pl, Zl) = /C^(p2> Z2\pl, Zl) 

22 

+ J dz J dplCv{p2, Z2\p, Z){-W{p, z))IC.o{p, Z\p^, Zl) + 

Zl 

Then after a simple algebra one can represent p3| ) in the form 



abs 



where 



BH 



dP^ _dP^^dP^ 

dx dx dx 



oo 



(46) 



dx 



-T -Re J dp J d^i J rf6^7(ei,6,a;)/C,(0,6|p,0) 



-co 

xcr{px)}Cy{p, 0|0, ^i) exp 



(47) 



dP° 



21 OO 



dx 



^Re J dzin{zi) J dz2n{z2) J dp^dp^ J d^i J d^2g{^i,^2,x)ICy{0,^2\p2, Z2) 

^(6-6 " 



-00 22 



(4J 



X(7{p2x)IC{p2, Z2\Pi, Zi)cr{pix)}Cy{pi, zi\0, ^1) exp 

Here T = dzn{z) is the optical thickness of the target (we assume that n{z) = at 
z < and z > L). The integrals over ^1^2 in (^7]), ( ^8)) of the products of the vacuum 
Green's functions and exponential phase factors can be expressed through the light-cone 
wave function \E'(x, p, Ag, K', A^) for the transition e e'7. At Ag' = Ag it is 



^(x,p, Ae, Ag/, A^) = 



1 I (y I 3 3 

[A^(2 - x) + 2Agx] ( — - iA^— ) j d^}Cy{p,0\0,^) 



X exp 



1 lax 



[A^(2 — x) + 2Aex] exp(— 'iA^v^)mgi^'i(pmea;) , (49) 



for Ag' = — Ag the only nonzero component is the one with A-y = 2Ae 

■ 



\E'(x,p, Ae,-Ae,2Ag 



V2 



ax-^ 



2/x, 



e'7 



d^}Cv{p, 0|0, ^) exp ( 7^ ) = — V 2ax^me_K'o {puiex) .(50) 
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Here Kq and Ki are the Bessel functions. Eqs. (^9]), (^) can be obtained by calculating 
the matrix element for transition e e'7 in time-ordered PT using the representation 
(§) for electron and photon wave functions. 

Making use of Eqs. (^9]), (^Oj) one can rewrite (^71) , (^Sf ) in the form 



/dp|^(x,p,{A.})|V(px), (51) 



1 

-^Re^ j dzin{zi) j dz2n{z2) J dp'^*{x, p, 
{Ai} 21 

X(T(px)$(x, p, {Ai}, 2:1, 2:2) exp 



21 

2(2:2 - ^i) 



(52) 



where 



^(x, p, {Ai}, zi, za) = y" rfp'/C(p, 22|p', ^i)^(a;, p', {A^}) a(p'a;) (53) 



is the solution of the Schrodinger equation with the boundary condition 
$(x, p, {Ai}, zi, zi) = ^(x, p, {Xi})a{px) . 

In Ref. p6| it was shown that the p_L-integrated cross section for a radiation process 
a ^ be can be written as 



da{a — > cb) 



j dpWl\x, p)a-abc{p) , (54) 



dx 

where W^'^{x,p) = |^a(a;, p)p is the light-cone probability distribution for transition 
a — > be, aabc is the total cross section of interaction with the target of dbe system. For 
the transition e —>■ e'7 the corresponding three-body cross section equals a{px). Con- 
sequently, the first term in ( ^61) equals the Bethe-Heitler cross section times the target 
optical thickness, i.e. it corresponds to the impulse approximation, while the second 
term describes LPM suppression. Thus, we have demonstrated that LPM suppression is 
equivalent to absorption for e"'"e~7 system. 

It is worth noting that at Lf ^ L the radiation rate for a composite target can also 
be represented in a form similar to that of Eq. (^4|). Indeed, in this limit the transverse 
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variable p is approximately frozen, and the Green's function can be written in the eikonal 
form 

a{pix) 



^(p2> ^2|Pi, zi) ^6{p2- Pi) exp 



dzn(z) 



Using (|55| ) we obtain in the frozen-size approximation from Eqs. (146|) , (|5lD , (|52D 

To-{px) 



= 2 J dp WPix, p)<l- exp 



(55) 



(56) 



Eq. ( ^6]) is analogous to the formula obtained in Ref. for the cross section of heavy 
quark production in hadron-nucleus collision. Within classical electrodynamics the LPM 
effect at ^ L was previously discussed in Ref. jTSl . 

Representation (^) has the virtue of bypassing calculation of the singular transverse 
Green's function. This renders it convenient for numerical calculations of the radiation 
rate for finite-size targets. 



2.4 Numerical results and comparison with the SLAC experi- 
ment 

For numerical calculations we need a more accurate parametrization of the dipole cross 
section, which takes into account the inelastic processes and the Coulomb correction. We 
write the dipole cross section in the form cr(p) = p^C(p) , where 

Cip) = Z^CM + ZCUp)- (57) 

Here the terms oc Z"^ and oc Z correspond to elastic and inelastic intermediate states in 
interaction of e~^e~ pair with an atom. Due to the steep decrease of the light-cone wave 
function \E'(x, p, {Aj}) at \p\ > l/rrieX the dominating values of p in (|5T|) are ~ l/rrie. For 
(^) they are even smaller due to the absorption effects. For this reason the probability 
of photon emission is only sensitive to the behavior of cr(p) at p < l/rrie ^ r^. In this 
region both the Cei and Cj„ can have only weak logarithmic dependence on p. This allows 



one to parametrize them in the form 

2 



Q{p) 



in 



P 



G/j \ Ct-i 



p 



Ana 



log 



2 a,; 



+ 



(1-27) 



(5^ 
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For elastic component Cei{p "^Ra) = CeK-^^)- We adjusted a^i and to reproduce the 
terms oc and oc Z, respectively, in the Bethe-Heitler cross section 



dx Smlx 



(4 - 4x + 3x^)[Z^(Fe^ - f{Za)) + ZF,^] + (1 - x) ^ ^ ^ | , (59) 
« log(184/Zi/=^) , ^ log(1194/Z=^/3) , 

f{y) = y'i: 



evaluated in the standard approach with realistic atomic formfactors This procedure 
gives a^i = 0.81 r^Z-i/^ exp(-/(Za)) and a^^ = 5.2rB^-2/3_ 

In Fig. 1 we compare the results of calculations (solid curve) of the bremsstrahlung 
rate with the one measured in for a gold target with L = 0.7%Xq ^ 0.023 mm 
(Xo is the radiation length) and 25 GeV electron beam. We also show the prediction of 
frozen-size approximation (^) (dashed curve), the radiation rate obtained for the infinite 
medium (long-dashed curve), and the Bethe-Heitler spectrum (dot-dashed curve). We 
have found that the normalization of the experimental spectrum disagrees a little with 
our theoretical prediction. The theoretical curves in Fig. 1 were multiplied by the factor 
1.03. This renormalization brings the calculated spectrum in very good agreement with 



the data of Ref. 22 



For 25 GeV electrons Lj ^ 0.94 ■ (lMeV/A;(MeV)) mm in the region of k shown in 

Fig. 1. One can conclude from this figure that the radiation density calculated using Eqs. 

(^), (PH), ( p2D is close to the prediction of the frozen-size approximation ( |56D for the 

photons with i^/ > 2L, while for the photons with < L it is close to the spectrum for 

the infinite medium. To illustrate the role of the finite target thickness better we present 

in Fig. 2 the LPM suppression factor defined as S" = ^pslfj"^^ as a function of the ratio 

h = L/Lf for several values of the photon momentum. The calculations were performed 
^ Recently we have analyzed the SLAC data including the multiphoton effects (iiep-pi/9805271). 
The results of this analysis are in very good agreement with the experimental data for all the targets used 
in ||2^. For the 0.7%Xo gold target the effect of multiphoton emission turns out to be small. It increases 
the normalization constant by ~ 3%. 
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for a gold target and 25 GeV electron beam. Fig. 2 demonstrates that the edge effects 
come into play at L <Lf. One can also see from Fig. 2 that for low photon momenta the 
edge effects vanish steeper. This fact is a consequence of a stronger suppression of the 
coherence length in radiation of soft photons. 

Fig. 2 shows that the suppression factor has a minimum a.t L ^ Lf for 100 and 400 
MeV photons. This minimum reflects the two-edge interference for a plate target. One can 
expect a more pronounced interference effects for structured targets. To illustrate the role 
of the interference effects in Fig. 3 we show our results for the LPM suppression factor for 
a two segment gold target. We have performed calculations for the same plate thicknesses 
and the gaps between two plates as in the recent paper by Blankenbecler [^]. The 
analysis ||19| was performed within the model proposed in Ref. [|1^], in which the medium 
was modelled by the potential U{p,z) = —p ■ E±{z) , where is a random transverse 



electric field. Qualitatively our results are similar to those of Ref. ||T9|. However, for 
our realistic electron-atom interaction the maxima and minima in the spectra are less 



pronounced than for the model medium used in Ref. |]19| . For a homogeneous target our 
spectrum differs from obtained by Blankenbecler by ~ 10 — 20%. 

The disagreement of our results with those of Blankenbecler is a consequence of im- 
possibility to simulate the Coulomb effects in the approach of Refs. ||18|, |19[. Indeed, 
using Eq. (|35D one can show that the model potential of Refs. [0, |19[ corresponds in our 
approach to the following choice of the dipole cross section 



jm^ f rfz(E^(0)-Ex(^)). 
n J 



Thus we see that in the approach of Ref. []T9| the Coulomb effects, leading to the important 
logarithmic p-dependence of the factor C{p) (|57D , are missed. We conclude that the model 
of Refs. |]rB|, is too crude for a quantitative simulation of the LPM effect in a real 
medium. 
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2.5 Probability of e^e pair production 

The probability of pair production by a high energy photon can be written in the form 
similar to Eq. (pSf). In this case the two-dimensional Hamiltonian reads 

H = ^ + v{p,i), (60) 

«(p.0 = -.^^i2^. (61) 

where /igg = E^x{l — x), x is the electron fractional light-cone momentum. The formation 
length for pair production is = 2E^x{l — x)/ml , and the vertex operator is given by 



9(^1, 6, x) = ^ ^ v(6) ■ v(6) + 77-/ , (62) 



where 

V(6) = Ve(^i) - Vg(^i) , 

Ve and Vg are the electron and positron transverse velocity operators. 

The light-cone wave function for transition 7 — > e~^e~ entering the representation 
analogous to Eq. (|^) is as follows: 

1 fa 

\E'(x, p, A^, Ae, Ag) = —\ - [A^(2x - 1) + 2Ae] exp{i\^ip)meKi{pme) , (63) 
/vr V z 

for Ae = — Ag, and the only nonzero component for Ag = Ag (in this case A^ = 2Ae) is 

^ 

^'(x, p, A^, Ae, Ag) = —V2ameKQ{pme) . (64) 
Ztt 

3 The LPM effect in QCD 

3.1 General expression for the probability of gluon emission 

Let us now consider the LPM effect for the induced gluon radiation from a fast quark. We 
discuss both cold nuclear matter and QGP. For QGP we use the GW model [|1^ treating 
QGP as a system of static scattering centers described by the Debye screened potential. 
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For the Debye color screening mass we use perturbative formula hd = (1 + np/GY^'^gs T 
28| , where Qs = ^^-nag is the QCD coupling constant, T is the temperature of QGP. 



Nucleons making up the cold nuclear matter are also treated as static scattering centers. 
Interaction of the fast quark and emitted gluon with each center will be described including 
one- and two-gluon exchanges. It should be noted that inclusion of the two-gluon exchange 
is absolutely necessary to ensure unitarity. 

The derivation of the gluon radiation rate follows closely the analysis of bremsstrahlung 
in QED. Similarly to Eq. the probability of gluon emission, Pg, is connected with the 
medium modification of the radiative correction to the probability to detect in the final 
state one quark, 5Pg, 

P, = -{SP, - SPr) ■ (65) 

Owing to the fact that —T* = Tq (here ^ are the color generators for a quark and 
an antiquark) the complex conjugated quark propagator is equivalent to the antiquark 
propagator. After summing over the final states of the target with the help of the closure 
relation 

l^f)(^fl = l, 

where is the wave function of the target after interaction with a fast quark, the 5Pq 
will involve only the diagonal matrix elements for the medium constituents. This means 
that only the diagrams involving color singlet (Pomeron) t-channel exchanges between the 
gg, qqg states and the medium constituents contribute to 5Pq. Consequently, in just the 
same way as in QED, we can obtain the expression for 5Pq in a medium introducing the 
corresponding absorption factor in the vacuum path integral formula for 5Pq. This allows 
one to obtain the formulas analogous to Eqs. (0), (|T^) and (|T8D . In the analogue of Eq. 
(|17|) the absorption factor contains the dipole cross section o"2 of interaction of qq pair with 
a medium constituent. In the absorption factor for the QCD analogue of Eq. ([T8| ) the 



corresponding cross section is the three-body cross section o"3 for qqg intermediate state. 
Namely this cross section enters the final formula for the radiation rate. For a quark 
incident on a target, it has a form that is analogous to equation ( |33D (we use notation 
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similar to that in the case of QED) 
dPg f f ^ ( iA^ 



dx 



2Re j c^ei/t^^exp -^ U(6,6,a;)[/C(0,6|0,ei)-/C,(0,6|0,6)] • (66) 



^(6, 6, X) = "-^^ v(6) ■ v(eo + . (67) 



Here the generahzation of the QED vertex operator (p^ to QCD reads 

as[4: - Ax + 2x^] 1a,m\x 
Tx "^^^)-"^^^) + ^ 

where 

and Vqi are the gluon and quark transverse velocity operators. The Hamiltonian for 
the Green's function /C is given by 

n=^ + vip,z), (68) 
n{z)(T3{p,x) 

v{p,z) = -t . (69) 

The Schrodinger masses are defined similarly to the case of photon radiation. The gluon 
formation length is 

2Egx{l - x) 
[m^a;^ + m^(l — x)\ 

rriq is the quark mass and rrig is the mass of radiated gluon. The latter plays the role of 
an infrared cutoff removing contribution of the long-wave gluon excitations which cannot 
propagate in the real nonperturbative QCD vacuum. In the case of QGP summation over 
triplet (quark) and octet (gluon) color states is implied on the right-hand side of Eq. (|69D. 
For a quark produced inside a medium through a hard mechanism the integration over 
^1 in Eq. ( |66D starts from the production point. Note that for gluon emission Lf ^ 
for soft {x 0) and hard {x — > 1) gluons. As a result, in both these limiting cases the 
Bethe-Heitler regime must obtain. 

In general case the three-body cross section for qqg state depends on the two transverse 
vectors: pgg and p^g, here = p^ — Pj. In terms of the dipole cross section it is given 

by m 



9 1 

'^^[Pqg^Pqg) = g [^^2 ( | P^g | ) + ( | P^-J )] - -a2{\pqq\) ■ (71) 
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However, in the case of interest tlie antiquark in tlie qqg system is located at tlie center- 
of-mass of the qg system, and 

Pqg = (1 - X)pqg , pqq = -XPgg . (72) 

For this reason the three-body cross section entering the imaginary potential (p9| ) can be 
written as 

9 1 

(^3ip,x) = -[o-2(p) + Cr2((l - X)p))] - -0-2(xp) , (73) 



where p = \pgg\ ■ Eqs. ([7^), ([73|) demonstrate that at x — * the color singlet qqg system 
interacts with medium constituents as octet-octet state, and as triplet-triplet state at 
X — s> 1. This is a direct consequence of the x-dependence of the transverse separations 
defined by Eq. ([7^). Notice that this makes evident that in the soft gluon limit one 
cannot neglect the transverse size of the qg system as was done in Refs. [|1T], pO[] . 
The dipole cross section can be written as 

a,{p)=C,{p)p\ (74) 

where C2{p) has a smooth (logarithmic) dependence on p at small p. For nucleon C2{p) 
in the small-p limit can be expressed through the gluon distribution 

C,{p) ^ f^xBg{xB, ^ 4) ' -4 ^ 10 ' (75) 
3 p'^ 

where x^ ~ {2xEqmpp^)~^ . For energies that are of interest from the practical view- 



point, the gluon density in Eq. (|75D can be estimated in the Born approximation, which 



corresponds to calculation of ct2 in the double gluon model of the Pomeron |3T[ . 

It is appropriate here to comment on gluon emission by a fast gluon. In this case 
qqg state will be replaced by ggg state. The ggg system can be in symmetric and anti- 
symmetric color states. As a result, the cross section for qqg state in the potential (|69|) 
will be replaced by the diffraction operator describing transitions between these two color 
states. However, in soft gluon limit the transition to symmetric color state can be ne- 
glected and we obtain the same Schrodinger equation as for gluon emission by a quark. 
The corresponding vertex operator is given by Eq. (^) times the color factor 9/4. 
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In this study we will use formula ( ]66| ) to evaluate the quark energy loss 



AE„ = E„ dxx 




(76) 



We will consider homogeneous nuclear matter and QGP. Of course, due to dependence of 
the probability of gluon emission on gluon and quark masses, our theoretical predictions 
are of approximate, estimating nature. Bearing this in mind, we will neglect the spin-flip 
transitions, which give a small contribution to the energy loss. Note that, in any case for 
a quark produced through a hard mechanism inclusion of spin effects in radiation without 
those at production vertex does not make sense. 

3.2 Gluon emission in an infinite medium in the oscillator ap- 
proximation 

Using Eq. (|54D for the transition q ^ qg one can show that the Bethe-Heitler cross 
section is dominated by the contribution from p < l/nig. In the case of gluon emission 
in a medium, the typical values of the transverse separations in the qqg system are still 
smaller due to absorption of the configurations with large transverse size . The smooth 
p-dependence of C2(p) at p < l/rrig allows one to evaluate the induced gluon radiation to 
a logarithmic accuracy replacing C2(p) by C2(pe//), where pe// is the typical size of the 
qqg system dominating the radiation rate (|6^). Then (T3(p, x) ~ C^{x)p'^, where 



Note that the large value of factor A in Eq. ([75|) is important from the viewpoint of 
applicability of the oscillator approximation. This allows one to use this approximation 
for a qualitative analysis of the induced gluon radiation even for a weak LPM effect when 

Peff ~ 1^9- 




(77) 



and the Hamiltonian ( |68D takes the oscillator form with the frequency 
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Using the oscillator Green's function p7|), we can obtain for the radiation rate per 
unit length 

dP da^^ 

-n^^S{7]), (78) 



dxdL dx 

where the suppression factor is defined by Eq. (PI]), and the Bethe-Heitler cross section 
is given by 

da'''' _ 4a,C3(a:)(4-4x + 2x^) 
dx 9TTx[mgX'^ + m^(l — x)] 

The dimensionless parameter t] in ([78|) reads 

(80) 

m^x"^ + mg[l — X) 

Note that the Bethe-Heitler cross section has the infrared ^ divergence. However, 
it is interesting that, in the limit of strong LPM suppression r] » 1, multiple scattering 
eliminates this divergence. Using the asymptotic formula ( ^3]) for 5'„,/(?7) at r/ >> 1, we 
can obtain from Eqs. (|78D , ([79|) in this regime 

dPg _ ag(4-4x + 2x^) 2nC3(x) 

(ixfiL ^ 37r ^ EgX^{l -x)' ^ ' 

The value of pe// in Eq- ( [TTD can be obtained from the diffusion relation pejj ~ 
{2A^eff /fJ'q'gY^'^- Here, as for the photon radiation, A^e// ~ = min(L/, 
This gives to a logarithmic accuracy pe// ~ [a;|n-E'gx(l — x)]^^/^. The elimination of 
the infrared divergence is a direct consequence of the medium modification of the gluon 
formation length. At rj >> 1 the medium- modified formation length L'j = Lf/rj « Lj, 
and the typical transverse size of virtual q'g system becomes small pe// << ^/^g- In this 
region the dynamics is scaling. As a result, the radiation rate (0) has only a logarithmic 
dependence on the gluon mass coming from the factor C3. Using the double gluon formula 
for the dipole cross cross section, we find from Eq. (|8TD at x << 1 for QGP 



dxdL ^ 37iEqX' 

where Ct is the second order Casimir invariant for the color center. For nuclear matter 
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after expressing C3 through gluon density, Eq. ( |8TD yields 



dxdL 



nasXBg{xB, 10/p' 



6Egx'^ 



13) 



Note that Eqs. ( ^2|) and (|83|) differ from predictions of Refs. and [|l7l by the factors 
V2/3 and 12, respectively. 

Ignoring the contributions to the energy loss from the two narrow regions near x ~ 
and X ~ 1, in which Eq. ( pT] ) is not valid, we find that, in the limit of strong LPM effect, 
the energy loss per unit length is 

dAE„ 



dL 



l.lasJnC3iO)E, 



q ■ 



4) 



3.3 Quark energy loss in hadron- nucleus collisions 

Let us now consider induced gluon radiation of a fast quark incident on a slab of nuclear 
matter of thickness L. This situation simulates gluon emission in hadron-nucleus colli- 
sions. From Eq. (|66D using the oscillator Green's function ( pTl) after some algebra the 
radiation rate can be represented in the form 



where 



d P da^^ 

L[m^x^ + m?g{l — x)] 



2Eqx{l - x) 



^5) 



6) 



I = L/Lj 

and rj is defined by Eq. (0). In terms of the dimensionless variables rj and /, the 
suppression factor S{ri, I) is given by 



2 111 yi 
S'^^\r],l) = T^Re J dyi J dy2exp 



i7) 



S^'^\r],l) = T^Re J dyi J dy2exp 


4> 



sin(0?/2) 

i{yi + y2. 



V 



cos(0?/i) (tan(0?/i) + 01/2 



^9) 
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S^^\r],l) = J dyi J dy2exp 



i{yi + 1/2 + Iri) 



1 



X 



(2/1 + y2 + ivY 



(90) 



Xl/i + 1/2) cos (0/?7) + (1 - (p'^ym) sin (0/7]) _ 
with = = exp(— z7r/4). The first term on the right-hand side of (|87D corresponds, 

in Eq. (|66D, to the contribution from the integration region ^1 < ^2 < L. The second term 
is associated with the region ^i<0<^2<-^, which gives the same contribution as the 
region < ^1 < L < ,^2- The last term in (|87D comes from the region .^i < and .^2 > L. 
The variables in (|g5D, (^), (0) in terms of those in (p^) are as follows: yi = (L — 
2/2 = (6 - in m, yi = {L- eOl^l, 2/2 = (6 - L)\n\ in (ID, = ~^,\n\, 

2/2 = (^2 — 1^1 in (0)- In deriving (^S]), (^9]), ( pO]) we have used a representation of 
the first Green's function in the square brackets in ( |66D in terms of a convolution of the 
oscillator and the vacuum Green's functions. At L ^ 00 the factors S*^^-* and S^^^ in 
Eq. (|87D vanish, while S^^^ tends to the infinite medium suppression factor (0). The 
finite-size effects come into play at L < L'f. For L « Lj Eqs. (|5), (^91), (|90D yield 
^(i'2)(r7,/) oc/2 and ^(3) ~ 1. 

In numerical calculations we take rrig = 0.75 GeV. This value of rrig was obtained in 



Ref. from the analysis of HERA data on structure function F2 within the dipole 
approach to the BFKL equation. It is also consistent with the nonperturbative 
estimate |Q of the gluon correlation radius in QCD vacuum Rc ~ 0.27 fm. Note that 
the hadronic size is bigger than 1/mg ^ Rc hj a. factor ~ 4 — 5. It is this circumstance 
that allows us to neglect the interference effects connected with gluon emission from 
different quarks. For real nuclei LPM suppression turns out to be relatively small. For 
this reason in Eq. ([77|) we take C2(Pe//) = C2(l/^g)- For scattering of the qqg system 
on a nucleon, we find from the double gluon model I^T] C2{l/mg) ~ 1.3 — 4 where the 
lower and upper bounds correspond to the t-channel gluon propagators with mass 0.75 
and 0.2 GeV, respectively. The latter choice allows one to reproduce the dipole cross 
section extracted from the data on vector meson electroproduction However, there 



is every indication |]32|, B3] that a considerable part of the dipole cross section obtained 
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in 



35| comes from the nonperturbative effects for which our approach is not justified. 



For this reason we take C2{l/mg) = 2, which seems to be a plausible estimate for the 



perturbative component of the dipole cross section ||3^. For quark mass, which controls 
the transverse size of the qqg system at a; ~ 1, we take rriq = 0.2 GeV. Notice that our 
predictions for AEg are insensitive to the value of m^. 

We performed calculations taking n = 0.15 fm~^ and ag = 1/2. Our numerical results 
in the region L < 10 fm can be parametrized in the form AEg ^ 0.1Eq{L/10{m)^ with 
/? ^ 0.9 - 1 for Eg < 50 GeV and p ^ 0.85 - 0.9 for Eg > 200 GeV. Our estimate is in 
a good agreement with the longitudinal energy flow measured in hard pA collisions with 



dijet flnal state |36| and the energy loss obtained from the analysis of the inclusive hadron 



spectra in hA interactions Note that our result differs drastically from the prediction 
by Brodsky and Hoyer ||]: AEg ^ 0.25(L/lfm) GeV. 

Our numerical calculations give the energy and L-dependence of AEg close to those for 
the Bethe-Heitler regime. This can be readily understood at a qualitative level. Indeed, 
for a quark incident on a target the radiation rate can be represented in the form analogous 
to Eq. (^) in QED. In the case of interest absorption effects at the longitudinal scale 
about the nucleus size play a marginal role due to small transverse size of the qqg system 
(~ l/rrig). As a result, the radiation rate must be close to the Bethe-Heitler one, and we 
immediately obtain AEg ~ EgLnasC^{Q)/m?g. Thus, we see that for real nucleus LPM 
suppression does not play an important role. Note that this clearly demonstrates that the 



approach that was used in Ref. [jT^ and which assumes strong LPM suppression is not 
applicable to hadron-nucleus collision. 



3.4 Energy loss of a quark produced inside a medium 

For a quark produced inside a medium the probability of gluon emission can also be 
written in the form (pSf). The suppression factor in this case is given by 

5(r/,0 = 5«(r^,0 + 5^'^(^,0, (91) 
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where S^^'"^^ are defined by Eqs. (]88|), (pQ]). From Eqs. (]88|), ( p9D one can obtain S'(?7, /) ~ 
— /^log/ for / << 1 . The physical mechanism behind this suppression of radiation at 
small L is obvious: the energetic quark produced through a hard mechanism loses the 
soft component of its gluon cloud and radiation at distances shorter than the time required 
for regeneration of the quark gluon field turns out to be suppressed. Notice that a similar 
suppression of photon radiation from an electron after a hard Coulomb scattering was 
discussed long ago by Feinberg |^ . 



Before presenting the numerical results, let us consider the energy loss at a qualitative 
level. We begin with the case of a sufficiently large Eg such that the maximum value 
of L'jt, Lj(max), is much bigger than L. Taking into account the finite-size suppression 
of radiation at Lj > L, we find that at high quark energy AE^ is dominated by the 
contribution from two narrow regions of x: 

x<5g^ Lmy2kEg , (1 - x) < 5, ^ Lml/2loEq , (92) 

where /q = min(l, I/77). In both the regions the finite-size effects are marginal and the 
energy loss can be estimated using the infinite medium suppression factor. For instance, 

AEgix <6,) ^ JdxS{r,{x),l = 00) . (93) 

3 

Using Eq. ( PDD one can show that ri{x <Sg)<latL< m^/2nC3(0). In this region of L 
in (|93|) we can put S{ri{x), I = 00) ^ 1 and find 

AEg ~ 0.25a,C3{0)nL'^ . (94) 



At L m^/2r?,C3(0) the typical values of rj in ( p3D are much bigger than unity, and using 
the asymptotic formula for the suppression factor we obtain 

AEg ~ a,C3(0)nL2 . (95) 

A similar analysis for x close to unity gives the contribution to AEg suppressed by the 
factor ~ 1/4 as compared to that for small x. Thus we see that at high energy AEg does 
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not depend on quark energy, and despite the infrared divergence of the Bethe- 

Heitler cross section has only a smooth m^-dependence originating from the factor C3. 
We emphasize that the above analysis of the origin of the leading contributions makes it 
evident that Li^ dependence of cannot be regarded as a direct consequence of LPM 
suppression of the radiation rate due to small angle multiple scattering. 

The finite-size effects can be neglected and AEq becomes proportional to L if 
Lj(max) <^ L. If in addition the typical values of r] are much bigger than unity, then the 
energy loss per unit length is given by formula (0). 

To study the infrared sensitivity of Ai?q, we performed numerical calculations for two 
values of mass of the radiated gluon = 0.75 and mg = 0.375 GeV. As in the case of 
a quark incident on a nucleus, we take rriq = 0.2 GeV and C2{peff) = C'2(l/'"^g)- ^^e 
case of QGP we take T = 250 MeV, and a = 1/3. For scattering of the qqg system on 
a quark and a gluon we use for C2(l/mg) predictions of the double gluon formula with 
the Debye screened gluon exchanges. In the region L < 5 fm our numerical results can 
be parametrized in the form 

The D and f3 as function of Eg are shown in Fig. 4 (nuclear matter) and Fig. 5 (QGP). 
In the region 5 < L < 10 fm /? in ( p6D is by 10-20 % smaller than for L < 5 fm. Note 
that Lj(max) ~ 5 — 10 fm for Eg ^ 10 — 40 GeV in the case of nuclear matter, and 
Eg ~ 150 — 600 GeV for QGP. Then from Figs. 4, 5 one can conclude that the onset of the 
1^2 YQgiin.e occurs at Lj(max)/L > 2. The closeness of f3 to unity at Eg ^ 10 GeV for QGP 
agrees with a small value of Lj(max) (~ 1 fm). Our results show that the m^-dependence 
of AEg becomes weak at Eg > 50 GeV. However, it is sizeable for ~ 10 — 20 GeV. 

Our predictions for AEg must be regarded as rough estimates with uncertainties of at 
least a factor of 2 in either direction. Nonetheless, rather large values of AEg obtained for 
QGP indicate that the jet quenching may be an important potential probe for formation of 
the deconfinement phase in AA collisions. A small quark energy loss obtained for nuclear 
matter indicates that the extraction of AEg from experimental data on deep inelastic 
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scattering on nuclei is a delicate problem. 



4 LPM suppression in hard reactions on nuclear tar- 
gets 

Another important example of the LPM effect in QCD is the well-known shadowing 
in hard reactions on nuclear targets. For instance, nuclear shadowing in deep inelastic 
scattering at small values of the Bjorken variable xb = Q"^ /2E^*mp, here and E^* are 
the photon virtuality and energy, respectively. This effect is similar to LPM suppression 
of pair production in QED. Calculation of the valence qq component of the shadowing 
correction Aa{'~f*A) = a{'~f*A) — Aa{'y*N) to 'y*A total cross section in the limit xb — > 



can be performed within the frozen-size approximation The light-cone path integral 
formalism allows one to take into account the parton transverse motion effects, which 
are important for evaluation of x^-dependence of nuclear shadowing. Nonetheless, an 
accurate analysis, requiring evaluation of medium effects for the higher qqgi...gn Fock 
states, is a difficult problem. However, within the Double-Leading-Log Approximation 
(DLLA) calculation of the leading twist contribution to Acr(7*y4) is greatly simplified. 

In the DLLA the parton light-cone variables and the transverse separations for the 
qqgi...gn state are ordered 

Xb << Xn « Xn-l « ■ ■ ■ « Xi << X < 1 , (97) 

^ « p2 « « . . . « p2 < _1 . (98) 

As a result, in calculating the leading twist shadowing correction, the subsystem 
qqgi-.-gn-i can be treated as a pointlike color-octet particle. Due to the ordering in 
the light-cone fractional momentum ( p7| ) the transverse motion of the center-of-mass of 
the qqgi...gn^i subsystem can be neglected, and only the motion of the softest gluon {g„) 
must be taken into account. Consequently, we can write Aa{'^*A) as 

A(j{^*A) = Aa^ai {l*A) + AasiP (7* A) , (99) 
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where Aa^a/ (7*^) corresponds to the qq Fock state of the virtual photon, while Aa^jp (7* A) 
gives the contribution associated with the higher qqgi...gn Fock states treated as a two- 
body octet-octet state. Both the terms on the right-hand side of (|99D can be written in 
the form similar to Eq. (^). For Aa^aiil*^) one can obtain (to simplify notation we do 
not indicate spin variables) 



Acr^ai(7*A) = —^^^^2 J dx j dh J dzin{h,zi) J dz2n{h, Z2) j dp'^1*,{x, p 

-00 Zl 

xa2(p)$^. (x, p, b, zi,Z2) exp 

where 



H^2 - Zl) 



(100) 



2E^*x{l - x) 



(101) 
(102) 



+ Q'^x{l — x) ' 
<l>^.(a;,p,b, 2:1,2:2) = J dp'fCq{p,Z2\p' ,zi)'^'i^*{x, p') cr2{p), 

n(h,z) is the nuclear density, is the photon virtuality, \E'^*(x, p) is the light-cone 
wave function for transition 7* qq. In Eq. ( p.02| ) JCg is the Green's function for the 
Hamiltonian 



where 



v{h,p,z) 



2fi 



v{h,p, z) 



.n{h,z)a2{\p\ 



(103) 



(104) 



and Hqq = E^*x{l — x). 

Using the light-cone wave functions for the qqgi...gn Fock states obtained in Ref. 
we can represent A(T3ip(7*y4) in the form 

Acr3ip(7*y4) = -- ^-^ Re J dxg j dh j dzin{h,zi) j dz2n{h,Z2) J dp X'^{xg,Q'^) 

i(z2 - Zl 



XB 
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x^^t {xg, p)o-2(p)$7* (a;^, p, b, 2:1, 2:2) exp 



(105) 



where 



jg _ 2Ey*Xg 

f ~ 2 — 

■' mi 
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^4x9,P,h,zi,Z2) = J dp'ICg{p, Z2\p' , zi)^*{xg, p') a2{p) . (106) 

The Hamiltonian for the Green's function ]Cg{p, Z2\p' , zi) in Eq. ( |106D can be obtained 
from Eqs. ( |103| ), ( p,04| ) replacing /igg by figgg = E^Xg, and a2 by |(T2. The factor (9/4)^ 
in Eq. ( |105D reflects the fact that the dipole cross section for octet-octet state equals 
(9/A)a2{p)- The factor in Eq. ( |105|) coming from the internal qqgi . . ■ Qn-i states is 
given by 

4/o(2v^) 



X\xg,Q') = ^I^E JdxJ dpp'a,{pWAx,p)\\ (107) 



where 



12, /a.("^r 



is the expansion parameter of the DLL A, and Io{z) ^ exp(2;) / is the Bessel function. 
The light-cone wave function describing the softest gluon entering Eqs. ( |105| ), ( |106| ) is 
given by p9| 



^^UXg.p) 



nig 



Ty/Xg 



G p G p 

Ki{mgpi)- — i - Ki{mgP2) ^ 



iPil IP2I 

where e is the gluon polarization vector, and Pi2 = p± r/2, |r| ~ l/<5- Note that, 
according to the derivation of Eqs. (|100|) , ( |105| ), the perturbative component of the dipole 
cross section entering these equations, must be evaluated in the Born approximation. 

Similar expressions can be obtained for shadowing corrections in Drell-Yan pair and 
heavy quark production. The results of numerical calculations of nuclear shadowing in 
hard reactions will be presented elsewhere. 



5 Conclusion 

We have discussed a new approach to the LPM effect in QED and QCD. This approach 
is based on the path integral representation of the light-cone wave functions. Using 
the unitarity we express the cross section of the radiation process a ^ 6c in terms of 
the radiative correction to the transverse propagator of particle a. Evaluation of the 
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cross section of transition a ^ be is reduced to solving the two-dimensional Schrodinger 
equation with an imaginary potential proportional to the total cross section of interaction 
of dbc state with a medium constituent. We have demonstrated a close relationship 
between LPM suppression for the radiation process a be and the absorption correction 
for abe state. 

For bremsstrahlung in QED we have evaluated the LPM effect for finite-size homoge- 
neous and structured targets. For structured targets we predict minima and maxima in 
the photon spectra. We have given a rigorous treatment of the Coulomb effects, which 
were previously treated only to a logarithmic accuracy. We have also included the inelas- 
tic process neglected in previous works. For the first time we have performed a rigorous 



theoretical analysis of the experimental data on the LPM effect obtained at SLAC p2 |. 



The theoretical predictions are in very good agreement with the spectrum measured at 
SLAC [^] for the homogeneous gold target with L = 0.7%Xo for 25 GeV electron beam. 

For the first time we have performed a rigorous analysis of the induced gluon radiation 
in cold nuclear matter and in QGP within GW model For a quark incident on a 



nucleus we predict AEg ^ 0.1£'g(L/10 fm)^, with (3 close to unity. For a sufficiently 
energetic quark produced inside a medium we find the radiative energy loss AEg oc L^, 
where L is the distance passed by the quark in the medium. It has a weak dependence on 
the initial quark energy. The dependence turns to as the quark energy decreases. 

We have also demonstrated that the developed theory of the LPM effect can be used 
for an accurate evaluation of the leading twist contribution to nuclear shadowing in hard 
reactions on heavy nuclei. 
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Figures 




k, MeV 



Figure 1: The bremsstrahlung spectrum for 25 GeV electrons incident on a gold target 
with a thickness of 0.7%Xo. The experimental data are from Ref. [22] The full curve 
shows our results obtained using Eqs. (46), (51), (52). The dashed curve was obtained 
in the frozen-size approximation (56) The long-dashed curve shows the spectrum for the 
infinite medium. The Bethe-Heitler spectrum is shown by the dot-dashed curve. 
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Figure 2: The LPM suppression factor for 25 GeV electron incident on a homogeneous 
gold target as a function of the ratio h = L/Lf and the photon momentum. 
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Figure 3: The LPM suppression factor for 25 GeV electron incident on a two segment 
gold target. The thickness of each segment is 0.35%Xo. The set of gaps is as follows: 
(sohd curve), 0.7%Xq (dotted curve), lA%Xo (dashed curve), 2.1%Xo (long-dashed 
curve), 3.5%Xo (dot-dashed curve). 
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Figure 4: The parameters D (a) and (3 (b) for the parametrization (96) for nuclear 
matter. The sohd hnes correspond to nig = 0.75 GeV, and the dashed ones to rrig = 0.375 
GeV. 
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